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DARK MATTER HALO MASS PROFILES 
Dan Coe^ 

Draft 

ABSTRACT 

I provide notes on the NEW, Einasto, Sersic, and other mass profiles which provide good fits to 
simulated dark matter halos (§31). I summarize various published c{AI) relations: halo concentration 
as a function of mass (SfTl). The definition of the virial radius is discussed and relations are given to 
convert Cmr, Myir, and r^ir between various defined values of the halo overdensity 

Subject headings: cosmology: dark matter 
— gravitational lensing 



galaxies: clusters: general — methods: data analysis 



1. MASS-CONCENTRATION RELATIONS 

The mass profiles of galaxy clusters appear to be 
more centrally concentrated than reahzed in simulations 



Cn (Broadhurst et al. 2008 Oguri et al. 2009 Sereno et al. 



20101) . It true , this may be evidence l or Early Dark En- 



ergy (see e.g., Grossi fc Springer2009 l) . Or per haps there 
is a less exciting explan ation (e.g.. Ba rkana fc Loeb 2009; 
Lapi & Cavaliere 2009 Meneghettfet al. 2010). For more 



details, see my discussion m Coe et al 



More 



conclusive results are expected from the CLASH HST 
MCT pro.iecf] and perhaps LoCuSS ( |Okabe et al.]|2009 



Richard et al. [2009 ) 



More massive halos generally have lower concentra- 
tions than less massive halos. This is seen in both simu- 
lations and observations, though less clearly so far in the 
latter (see below). More massive halos form later, re- 
sulting in lower concentrations refiecting the lower back- 
ground density at the time of formation [Navarro et al.| 
pl96j). 

Lot a given radial mass profile (see qSl), the concentra- 
tion is defined as: 



'^vir I '^- 



(1) 



a mishmash ratio of the virial radius r.yjr and the radius 
r_2 at which p ex r~^ . Eor an NEW profile, r_2 = r^. 
The definition of the virial radius ryir is discussed at 
length in ^ but it is typically approximated as the re- 
gion within which there is an average overdensity of a 
certain value (Ac ^ 100 or 200) above Pcrit- For clar- 
ity, one may quote the exact value of A^, used: C200, for 
example. 

In principle concentrations could be derived using any 
radial fitting profile (Q. However the choice does mat- 
ter as the profiles behave differently between r_2 and 
Tyir- Concentrations de rived from NFW and Einasto fits 

are compared in 



2008) 



to the same halos (Duffy et al. 

Fig. [SJ Einasto c{M) relations have been d erived for the 
Millennium simulation r elaxed ( Gao et al.„2008} and all 
Hayashi fc White l ( [2008 [) halos. Below we locus on c{M) 
fr 



relations derived from NEW fits 

^ Jet Propulsion Laboratory, California Institute of Technol- 
ogy, 4800 Oak Grove Dr, MS 169-327, Pasadena, CA 91109 
^ Cluster Lensing And Supernova survey with Hubble, 
|http://www.stsci.edu/~postman/CLASH7^ 



1.1. Current c{M) measurements from NFW profile fits 
The current best estimates for c {M, z) are probably 



those given by Duffy et al. ( 2008 ) and [Maccio et al. 
(20081. Their findings are similar. Both analyze sim- 
ulations which use the WMAP5 cosmology, resulting in 
~ 20% lower concentrations than W MAPl (Tabl e [Sl as 



used in the Millennium simulation (Neto et al 




for example. Duffy et al. (2008) find that present-day 
(z — 0) halos follow the following mass-concentration re- 
lation: 



C200 =i 5.74 



M, 



200 



2 X lO^'^h-^Mp 



-0.097 



(2) 



They provide a separate relation for relaxed clusters 
which are more symmetric and thus better fit by radial 
profiles such as NFW. These have 15 — 20% higher con- 
centrations (Fig. [7]): 



C200 - 6.67 



A/2, 



2 X lO^'^h-^Mp 



-0.092 



(3) 



Intrinsic scatters are A log 2^0(^200) — 0.15. These rela- 



tions are pl otted in Fig. ^ along with corresponding re- 
la tions from Maccio et a^ (j2008 ) . 

Duffy et al. ( |2008 ) also supply fitted functions to halos 
spanning the redshift range z = — 2. Full: 



C200 



and relaxed: 



5.71 



(1 + ^) 



6.71 



0.47 



M2, 



2 X 1012/i-lMfr 



-0.097 



(4) 



C200 



M- 



200 



(l + z)"-44 \2x 1012/i-iMe 



-0.092 



(5) 



In their Table 1, they provide uncertainties for these fit 
parameters as well as corresponding values for Cmr and 
Myir- These C2oo(-A/200: z) relations are plotted in Figs 
& |4J Also plotted ar e the correspondin g Cvir{Myir,z^ 
relations provided by (Duffy et al. 2008). In Fig. Ill we 
plot the Bullock et al. (2001) c ex JT+z)~^ scaling for 
comparison. Note that Dutty08 find weaker dependencies 



on redshift: C200 oc (1 



-0.45. 



ex (1 + z)- 



-0.70 



Halo concentrations are sensitive to cosmology. A 
higher a% causes halos to form earlier, resulting in higher 
concentrations. This was the case in the Millennium 
simulations which used the WMAP 1-year cosmology, 
including erg — 0-9- This yields concentrations ^ 20% 
higher than found in simulation s which use WMAP5's 
o-g = 0.79^ ( [Duffy et al.]|2008| ). The effect of cosmol - 
ogy was explored in more detail by [Maccio et al. (2008). 
These effects are shown in Fig. [5] 

Various derived c{M) relations (for z = 0) are plot- 
ted in Fig. |6| In our Tables [31 [4l and [5l we provide , . 
C20o(M20o), c„,AM,„ ), and cTlCz) re la tions, respec- weighmg m at M ~ 10 M^/h 



But the expected concentration is actually lower (c2oo ~ 
3.0, exacerbating the disagreement with observations) for 
four reasons: 

• C200 < Cvir 

• z > 

• WMAP7 vs. WMAPl 



M^ir sa 1.4 X W^Mq/H > W^Mq 



We also note that Bullock et al.| ([2001 1 did not simulate 
halos as massive as A1589, with their rnost massive halos 



tively^ as derived b y Duffy et al. (120 081), iMaccio et al 
([20081, [Neto et al.[ ^7\, [BuUocre t al. (2001),~ [HeS'- 
nawi et aH ( [2007[), arid fCentile et al.[ ( |2007 )[^ The latter 
is the original JNli'W c{M) prescription updated to the 
WMAP3 cosmology. 

The various simulations considered here are outlined in 
Table^ The relevant details of their adopted cosmolo- 
gies (SXri, o-g) are given in Table [8[ We provide the range 
of halo masses produced in eacn simulation. The dan- 
gers of extrapolating c{M) relations beyond these ranges 
h ave been cit e d by^Zhao et al._(2003J , for example. 

'Neto et al.[ ([20071 their I'ig. "syimd that 10,000 parti 



cles within the virial radius are required to yield robust 
concentration measurements. They note that using fewer 
particles introduces scatter but does not appear to intro- 
duce bias in their concentration measurements. 



jHennawi et al. ( 2007 1 measure significantly larger con- 
centrations tor galaxy clusters in their simulations. Their 
cluster concentrations are ~ 50% and ^ 80% larger than 



found by Duffy et_ al 
respectively (i'ig. pf 



([2008^ 



right) . 



and Maccio et al 
rhcir use ot erg 



(2008), 
USS^rob- 

20% 



ably only results m concentrations inflated by 
compared to the WMAP5 cg = 0.796 simulations. The 
remaining disagreement may be a result of their halo den- 
sity fitting procedure which they claim is better for com- 
parison with lensing measurements. Specifically, they as- 
sign large uncertainties to radial bins with large subhalos. 
This may bias the fitted profiles to be low at large ra- 
dius (where large subhalos typically reside) resulting in 
higher concentrations. These results may considerably 
ease tensions between observed and simulated halo con- 
centrations. The differences in fitting procedure should 
be better studied and understood. 

1.2. Care in citing concentration expectations 

A concern often noted is that the concentration mea- 
sured for A1689 (in every study to date) is higher than 
that found in simulations for a halo of A1689's mass. The 
concentration found in simulations has bee n cited loosely 
asc~5orc~5.5 using a relation given by [Bullocket al.[ 



(2001) 



9 



M„ 



-0.13 



(6) 



1 + z Vl-3 X lOi3/i-iAf0^ 

^ This value is in excellent agreeme nt with the WMAP 7-vear 
maximum likelihood value erg = 0.803 (iKomatsu et al.|2010||. 

* I'll take a crack at completeness, and mention other pa pers 
with cjM'] relations neglected here (for no particular reason): |Eke | 
let al. (2001); Wechsler et al. (2002); Alam ct al. (2002); Zhao et al. 
fe03); Uolag ct al. (2004); Kuhlcn ct al. (2U05); Lu ct al. (2000); 
[SFaw et al.| ( |2006p ; |Dutton et ai.| l |2U07p ; |Gnedm et al.| ( |2007p . And 
1 nave probably missed still otnersl 



1.3. The observed c[M) relation 

Based on a compilation of 62 published measurements 
(including 10 new measuremen ts) of halo Cmr and M^r, 
Comerford & Natarajan ( 2007 ) find the following relation 
(with a large scatter): 



14.5 ±6.4 

l + z 



M„ 



1.3 X 1013/i-iM 







-0.15±0.13 



(7) 



For clusters as massive as A1689, Comerford' s relation 
converges toward that of Hennawi et al. ( 2007 ) , the for- 
mer being only slightly higher. 

This and other observed c{M) relations are shown in 
Fig , ^and detailed in Table[6J The [Comerford fc Natara- 
jan (2007) compilation includes both lensing and X-ray 
determinatio ns of c and M , inclu ding the X-ray samples 
presen ted by Buote et al.j ( 2007 ) and Schmidt & Allen 
(20071. Each ot these papers presented their own c{M) 
relation. A recent c{M) relation from weak lensing o f 



individual halos was presented by Okabe et al. (12009') 



And c{M) derived from stacked w e ak le nsing analyses 



5 

th, 



al. (2007| and Mandel 



were presented by Johnston et 

baum et al. (2008). It seems apparent that one should 
study a wide enough range of halo masses to obtain a 
confident c{M) relation. 



2. OVERDENSITY WITHIN THE VIRIAL RADIUS 

Various conventions are used to define the virial mass 
and radius. We explain and show how to convert between 
different definitions. 

2.1. Overdensity Definitions 

The virial radius r^r designates the edge of the halo. 
Within this radius, objects are supposed to be "virial- 
ized" : gravitationally bound and settled into regular or- 
bits. Outside this radius, objects are not in orbit al- 
though they may still be infalling. In practice, there is 
no sharp dividing line the two regions. And even if there 
were, it would be extremely difficult to discern observa- 
tionally for a given massive body. Meanwhile, the ob- 
jects we study are not always virialized. In fact, galaxy 
clusters are the largest bodies which have had time to 
virialize given the age of the universe. Thus some of 
the clusters we observe have virialized just recently, but 
many are still in the process of doing so. 

Despite these complications, we can define a virial ra- 
dius for a massive body based on theory and simula- 
tions. Early theoretical work ( jPeebleS|il980| ) predicted 



that a sphere of material will collapse if its density ex- 
ceeds 1.686(1 + z) times that of the background. After it 
collapses and virializes, the sphere will obtain an average 
density 

Ac « IStt^ w 178 (8) 

times the critical density Pcrit{z) at that redshift, where 



Pcrit 



3H^z) 
8nG 



(9) 



Cole fc Lacey ( 1996 1 cited this as a theoretical result and 
then confirmed it in simulationsp 



of^ 



Navarro et al. ( 1996 1 adopted the nice round number 
200, which has been used commonly ever since 



to allow for easy comparison between papers. But the 
Ac ~ 178 result was obtained in an Einstein de-Sitter 
cosmology of {flm, ^a) = (1, 0). In the concordance 
cosmology (Jim, ^a) — (0.3, 0.7), we find a much lower 
value of Ac ~ 100, as we describe next. 

At least three different forms have been given for Ac as 
a f unction of cosmology. Fo r a flat universe {flm + ^A = 
l),|Bryan & Norman|(1998| give 



IStt^ - 82r2A 



39f}i. 



(10) 



An approximation to this is given as (Eke et aL]|1998) 



Ac w usa 



0.45 



(11) 



And Nakamura fc Suto ( 1997 their Eq. C19; see also 
IfaIFy||20U0[ their Eg. AlTJgwe 



Ac « 187r2(l + 0.4093a;2-^i^2)f7,; 



(12) 



with X = {l/^m.o - 1)^/^(1 + zy^ and ilmiz) = 
1/(1 -f- x^). Given the current concordance cosmol- 
ogy with Qrn = 0.3, these different expressions yield 
Ac — 101.1, 103.5, 100.3, respectively for a halo at z = 0. 
Or given 17„ = 0.25, Ac = 94.2,95.4,93.5. We note 
^m = 0.25 is in better agreeme nt with the WMAP 7- 
year value ([Kom atsu et al.|[2010[ ) and h = 0.742 ± 0.036 



from Riess et al 



( .2009^ . 



The overdensity is often quoted as a factor A^ 
the mean background density pm — ^mPcrit'- 



Phalo — AcPc 



A. 



Pr, 



virHin 



above 



(13) 



et al. 



et al 



With Ac = Ayir^jn, Ac = 101.1 corresponds to A„i^ = 
337 fo r r^m = 0. 3. This value is cited by e.g., [Bullock 
20011 ) and 'Grah am et all (|2006[). Using the 'Eke 
1998p relation and the Spergei etaL ( 2003 ) first- 
WKP value of Q„, = 0.268, Merritt et al.l (120061 



year 

among others) give A^ij, — 368. 

To facilitate comparison among current and future 
investigations, we propose that a value of Ac = 100 
be adopted for present-day halos. This corresponds to 
Ayir = 333 given flm — 0.3, or the nice round num- 



for fim = 
yields Ac 



400 gi ven fj^ 
25, the 
94 A 



0.25. We also note that 



0.25, the Bryan & Norman (19981 expression 
376. 



^ [Cole fc Laceyl | |1996| spoke of this factor as the overdensity 
above pm = ^imperii rather than above pcrit- But as Qrn = 1 
in their simulations, the twro densities were equal and thus inter- 
changeable. 



While results from simulations are most often reported 
for present-day halos. Nature provides us observers with 
images of clusters as they were in the past. Thus in 
the expressions above, we should replace the present day 
values of fi^.o and ftA,o (here "0" subscripts have been 
added for clarity) with: 



i IjjiyZ ) 



1 



l + (r!A,0Mn,0)(l + ;2)-3 



(14) 



and flA{z) = 1 — rim(z). For the massive galaxy cluster 
A 1689 aX z = 0.1862 and adopting Um = 0.3, the widely 



used Bryan fc Normaii] (|1998[) expres sion yields Ac — 
116.6 and theliNakamura & Suto|(|1997|) exr 



( 1997 ) expression yields 

Ac = 115 . The latter was adopted by Broadhurst et al. 
( 2005a| b private communication) so we adopt it as well 
for consistency in C oe e t al. (2010). 
In Figs. 13 14 and |l5[ we plot ri„i(z), Ac(2;), and 

^vir \Z ) . 

2.2. Conversion between overdensity values Ac 

If the mass profile is well described by an NEW pro- 
file, then it is straightforward to convert Cyir^ i^vir, and 
Myir between different conventions of Ac (c.f.. Fig. Isl). 
Converting from C200 (Ac = 200) to cus (Ac = 115), lor 
example, simply involves findin g th at value of C115 which 
yields the same value of 6c (Eq. 38) as did C2oo- This can 



be accomplished by a simple root finding program, but 
the relation is very linear as shown in Fig. [91 Here we 
provide expressions 

C94«1.328c2oo+ 0.272 (15) 

cioo«1.298c2oo+ 0.246 (16) 

cii5«1.232c2oo+ 0.189 (17) 

which are accurate to within 0.5% for 2 < C200 < 25. 
These factors may be generalized: 

Cyir w ac2oo + b (18) 



aft 


-1.1191ogio Ac -H 3.537 


(19) 


6ft 


-0.9671ogio Ac -H 2.181 


(20) 



to yield Cyir within 1% for 3 < C200 and 70 < Ac < 140. 
These factors are plotted in Fig. [TO] 

An alternate expression gives Cyir as a function of C200 
and Ac: 



C200 



„0.9 
^200 



10^ 



(21) 



-5^ A 1.82 



p = -(8.683 X 10"^)A, 



(22) 

to within 1% for 3 < C200 < 35 and 85 < Ac < 165 (z < 1 
or so). 



See also Hu & Kravtsov (2003 Appendix C). 



2.3. Virial Mass 
Virial mass (the mass within ryir) is given by 



My 



4 .3 . „ . ^ rl^AcH^z) 

-7Tryi^AcPcrtt{z) = ^ 



(23) 



M„ 



In Fig. 12 we plot this simple relation between ryir and 



Quoted values for M200 and r2oo can also be converted 
to Myjr and ryir as a function of C200 and Ac as plotted 
For a given NFW curve (with fixed Vg and 
simply scales with Cvir since r^ stays fixed. So 
— c„ir/c2oo. As for the M„jr conversions, we 



in Fig 

r„ir/ 1^200 



11 



solved for those numerically. 

3. MASS PROFILES 
3.1. Double Power Laws 
In dark matter simulati ons, galaxy and cluster halos 



Navarro et al. (1996 1997) were all shown to have mass 



density profiles well approximated by the NFW profile: 



p{r) 



(r/r,)(l + r/r,)2 



(24) 



This profile behaves as p ex r^^ in the core, p ex r~^ at 
r = r^, and steepens to p ex r~^ in the outskirts. The two 
fit parameters ps and r^ were shown to be related and a 
function of halo mass. This "universal" profile is still a 
good approximation to today's simulated halos. However 
the higher resolution does reveal subtle differences. 

Deviations were sought for using a generalized version 
of the NFW profile (Hernquist 1990, his Eq. 43; see also 
Zh^p96| IWyithe et al.||lJ001| ): 



p{r) 



2(/^-7)/" p^ 



(r/r,)''[l + (r/r^)"]('3-7)A 



This profile behaves as p ex r '''in the core, and p ex 

in the outskirts. The rate of transition is g overned by a. 

Where NFW found ' ' ' 

instead 



(1999) 



(a, /3, 7) 
found best fits 



(1, 3, 1), [ Moore et al 
_ of (1.5 

tantly, the inner profile was steeper: 7 



3, 1.5). Impor- 



1.5 



= 1.5, p ex r" 

There were many other attempts to accurately resolve 
and measure this inner slope, including [Diemand et al. 
(p005) who found p (x r~ ^ ^ 



-1.2 



^rhe fully generalized form in Equation 25 proves a bit 
too gener al with large degene racies between the free pa- 
rameters (Klypin et al. 2001). Thus, in their efforts to 
determine the central slope 7, authors often use one of 
two constrained versions of Equation [25j either a "gen- 
eralized NFW" profile with (a, /3, 7) ^1, 3, 7): 



p{r) 



23-7 



Ps 



(r/r3)7[l + (r/r,)](3-7)' 



(26) 



or what we might call a "generalized Moore" profilqjwith 

(a, /3, 7) = (3-7, 3, 7): 



p{r) 



2ps 



(r/r,)'>'[l + (r/r,)3-7]- 



(27) 



Meanwhile, Dehnen & McLaughlin (2005) found 

(a, /3, 7) = U/i), 31/9, 7/9): 



p{r) 



2pl 



(r/r,)V9[l + (r/r,)4/ 



916 



(28) 



° This latter form is also often referred to as a "generalized 
NFW" profile, although strictly speaking it can only exactly re- 
produce the Moore profile and not that of NFW. 



and when accounting for anisotropy, the more general 

{a, /?, 7} = {(3-7)/5, (18-7)/5, 7}: 



p{r) 



2pf 

(r/r,)7[l-H(r/r,)(3-T)/5]6- 



(29) 



3.2. Continuously Varying Power Laws 



The original NFW proponents proposed a new profile 
which gradu ally flattens all the way toward t he center 
([Navarro et al. 2004 K This profile was found (Navarro 



eraT:^004; Merritt et al. 2005 2006) to yield better 



fits to a wide range ot simulated dark ma tter halos than 
did the generalized NFW profile (Eq. [26]), which has an 
equal number (3) of free parameters, including the cen 
tral slope. Inner slopes as steep as p{r) '^ ---i^ 



are 



clearl y ruled out by recent simulations (Navarro et al. 
20101) 



The new Navarro et al.[(|2004 ) fittin g form \ yas qu ickly 



recognized ( [Merritt et al.|20U 5') as the lSersic ( 1968) pro- 
file generally applied to fitting the light distributions of 
elliptical galaxies. The implications are intriguing: that 
the collapse of massive bodies, be they luminous or dark 

matter, may lead to sim ilar profiles. 

However to be precise, Navarro et al. ( 2004 1 fit a Sersic- 
like profile to 3-D density distributions, where the Sersic 
profile was fit to 2-D surface density distributions (of 



(25) light). Einasto (1965) was first to use such a density 



law to describe a 3-D distribution, namely the spatial 
distribution of old stars within the Milky Way. 

Today we distinguish between the "Einasto" and 
"Sersic" mass profiles. The former is fit to 3-D mass 
density p{r) while the latter is fit to 2-D projected mass 
distributions S(i?). Projected and deprojected approxi- 
mations to the Einasto and Sersicprofiles, respectively, 
have a lso been deri ved (see Table [l|. 

The Sersic (1968) profile is given by: 



nm 



, exp' 



l/ri 



(30) 



There are three free parameters: Se, Re, and n, with 6„ 
being a function of n (given below) such that half the 
mass is contained within R^. Note that the total mass of 
a Sersic profile is finite, unlike that for an N FW profile. 
A 3-D deprojected approximation is given by [Prugniel &[ 
Simien[(ll997[ ). 

'I'he Einasto mass profile is a similar function but of 
3-D mass density p{r): 



p{r) 



p_2 exp 



r 
r-2 



1 



(31) 



where p_2 and r_2 are the density and radius at which 
"^. The concentration is defined as Cmr = 

0.17 for 



p(r) ex 

Tvir/'r-2- Navarro et al. (2010) found a 
galaxy -sized halos in the Aquarius simulation. |Gao et aL] 
( 2008 1 concur and found a increases to ^ 0.3 for th e most 
massive clusters in the Millennium simulation. [Duffy] 



et al. (2008) reduce the Einasto profile to two free pa- 
rameter by using the "peak height" a{v) relation from 



Gao et al. (2008). A 2-D p rojected appro ximation of 



Einasto is given by Dhar & Williams (2010) 



TABLE 1 
Sersic-like profiles of 2-D and 3-D density 



3-D p{r) 



2-D E(i?) 



Einasto 



(19651 



projected 



Dhar & Williams 



(20101 



Prugniel & Simien ( 1997 1 deprojected 



Sersic 1 1968 1 



3.4. Profile Details 

Here we provide useful expressions derived from the 
NFW and Sersic profiles. 

3.4.1. NFW Profile 

Simulated galaxy and cluster halos [Navarro et al.| 
(11996 1997) were shown to all have mass density pro- 
files well approximated by the NFW profile: 



TABLE 2 

3-PARAMETER FITS: DEVI ATIONS FROM HAL O 
PROFILES MEASURED IN IMeRRITT ET AlI 

( |2006[ THEIR Table 4) 



Model 


6 clusters 


4 galaxies 


Einasto 

Prugniel-Simien^ 
generalized NFW 
Dehnen-McLauglilin 


0.028 
0.025 
0.032 
0.034 


0.026 
0.030 
0.028 
0.023 


^ deprojected Sersic 







Merritt et al.| (2005J |2006[ ) experimented with both of 
these and other fits to 3-U mass density profiles of sim- 
ulated halos. In the latter paper, they compared the 
performance of various formulae fit to 10 simulated ha- 
los (6 cluster-sized and 4 galaxy-sized). Their results 
for 3-parameter fits are reprinted here in Table [21 The 
four tested profiles yielded similar results. Einasto per- 
formed a bit better across the board. Prugniel-Simien 
(the deprojected Sersic profile) performed a bit better 
for clusters. Dehnen-McLaughlin performed a bit better 
for galaxies. And the generalized NFW profile (Eq. 26 1 
was not far behind. 

Given the small number of halos tested, and the sim- 
ilarity of the performances, all of these fitting formulae 
might still be considered reasonable choices. However 
the Einasto profile has become e speciall y popular (e.g., 
Hayashi & White 2008; Gao et al . 2008; D uffy et al.|2008l 
TTavarro et al.||2010p . 



Recently a new fitting formula was proposed by 



et al. ( 2009 ) . It yields superior fits to the two high reso 



lution halos tested (Via Lactea 2 and GHALO): 



(32) 



with A ~ 0.1. 



3.3. Power Law Behaviors 



Power law slopes for the above fitting formulae as a 



function of r are plotted in Figs. 



fTe) [iTl [ll The slopes 
were calculated numerically as dmp/dlnr. For fun, we 
note this is equivalent to {dp/dr){r/p). The radii r are 
given in units of r_2 at which p ex r^^. 

We see that in principle, we should be able to dis- 
tinguish between these various profiles in both observed 
and simulated halos given sufficient resolution at a large 
enough range of radii. Such clear determinations have so 
far eluded us. 



In Fig. 19 we plot power law slopes for p(r), k{R), and 
M(< i?) Tot the NFW profile. 



p(r) 



{r/rs)il + r/r,y 



(33) 



The two fit parameters Ps and r^ were shown to be re- 
lated and a function of halo mass, as we discuss below. 
But the parameter making all the buzz these days is the 
central mass concentration: 



^vir / ^ 



(34) 



where ryir is the virial radius of the mass halo. As dis- 
cussed above, the virial radius is estimated as that which 
contains an average density AcPcrit, for a total virial 
mass of 






(35) 



For an NFW halo, the mass within a sphere with 
radius r = xr^ c an b e found by simply integrating the 
NFW profile (Eq.lsl: 



(36) 

(37) 



M(r)=47rrf / dx'x' 
Jo 



x'{l + x')^ 

4TTpsrU\n{l + x)--^ 
\ I + x 



Combining Eqs. 34 351 and ^71 we find that the concen- 



tration parameter c c an be obtained from t he following 
expression, as given in 



Navarro et al. (1996): 



Pcr^t ^" 3 In(l + C)-C/(1+C)' ^^^' 

To fit the NFW profile to our gravitational lensing 
mass maps which measure projected surface density, we 
integrate the NEW p rofile along the line of sight (e.g., 
Golse & Kneib 2002) to find the projected surface den- 
sSyl 



E(i?) = 2psrsFiX) 



(39) 



with R = Xr, and 



1 



1 



X2-1 



Vi-x- 



: cosh 



X 



{X<1) 

(40) 



F{X) = < 

.x^-iy^ VX2 - 1 """ x^ 

Integrating once more over the area within R, we find 



^ : cos-i - ) {X >l) 



6 



the total mass within a cylinder of radius R 

M{R) = 47rr3p,G(X) (41) 



with 



G{X)=\n^ + { 



, cosh"^ — (X<1) 

Vi-^2 X'' ' 



I Vx 



1 -1 1 

- cos — 
1 X 



{X = l) (42) 
(X>1) 



This should not be confused with Eq. [37] which gives mass 
within a sphere of radius r. 

From this we can obtain the shear due to an NFW 
mass profile: j{R) = k{R) — k{R): 



7(i?) = 2.. (^ - Fix) 



(43) 



The quantity measured in weak lensing studies is the 
reduced shear: 



9 



jDLs/Dsh 
1-{Dls/Ds)k' 



(44) 



where we have finally given the redshift dependence. 
(All previous expressions were given for a fiducial Icnsed 
source at Zs = oo.) 

3.4.2. Sersic Profile 



We now give the Sersic 



derived from it (e.g., |(jra 



(|1968[) profile and quantities 



bam & Driverir2UU5l iTerzic & 



Graham|2005 1 . Note that the Kersic profile is commonly 
used to describe the (projected 2-D) hght profiles of cl 
liptical galaxies. Here instead it will be discussed as de 
scribi ng projected m ass profiles. 
The Sersic ( 1968 1 profile is given by: 



i;(i?) == EeCxp ■! -bn 



l/n 



-1 



(45) 



There are three free parameters: Sg, i?e, and n, with 
bn being a function of n (given below) such that half the 
mass is contained within Rg. (Note that the total mass of 
a Sersic profile is finite, unlike that for an NFW profile.) 
The total projected mass within a radius R is given as: 



M{R) = 2TT^eRlne''^b-^''j(2n,x) 



where 



x^b„{R/R, 



xl/n 



(46) 



(47) 



and j{a,x) — L dte~*t°-~^ is the incomplete gamma 
function (with the "hat" used to distinguish 7 from the 
lensing shear 7). Thus to satisfy M{Re) = \M{R = 00), 
bn must obey: 



r(2?i) = 27(2n,6„) 



(48) 



where T(a) = 7(0, cx)) is the complete gamma function. 
In SciPy's "special" package, we find a routine to quickly 



calculate 6„ — p;an'mai ncinv(2 * n, 0.5). An app roxima- 
tion may also be used ( Prugniel fc Simien||1997 ): 

bn~2n- 1/3 + 0.009876/n (49) 

Lensing properties of the Sersic pr ofile have been de- 
rived and exp lored in Cardone (2004) and Eliasdottir & 
Moller (2007). Of special mterest here is the weak shear 
7 = K — K. The average k, within R can be derived 
straightforwardly from the above expression for M{R): 



k{R) 



M{R) 



2K.f,ne 



''"a;-2»^(2n,x 



T:R^Y,crit 

where we have introduced Kg = Y,e/T,crit 
k(R) = Y,{R) /Ticrit can be rewritten as: 

Thus we find 7(i?) = R{R) - k{R): 



(50) 
Meanwhile, 

(51) 



7(i?) = Kec''" {2nx-^"j{2n,x) - e"^) 



(52) 



with the reduced shear given as g = {'^Dls/ Ds)/[l — 
kDls/Ds). 

There are fewer published fits of S ersic profiles to simu- 
lated cluster halos. We do note that Merritt et al. (2005 ) 
found n = 2.38 ± 0.25 for their cluster sample. 



We thank Angelo Neto for useful conversations about 
the Millennium simulation and their study of halo pro- 
files. This work was carried out at Jet Propulsion Lab- 
oratory, California Institute of Technology, under a con- 
tract with NASA. 
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Fig. 1. — For fixed M200 = lO-^^M©, we show what an "over-concentrated" looks like: C200 = 9 versus the expected concentration 
(c200 = 3). This is roughly the case for A1689 (Coc ct al. 20101. Dots mark r^ and dashed lines mark r^i^ and Mvir- Thicker lines are 
used for two decades in radius about rs to emphasize the shifting of the curves. 
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Fig. 2. — Expected NFW concentration C200 ^'Hd l-c scatter as a function of halo mass M200 for a-U clusters (left) and relaxed clusters 
(right). 
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Fig. 3.— Expected NFW concentration c(M) for ^ = halos and c{M, z) for 2 = 0, 1, 2 halos from [Duffy et al.lpOOS} . 
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Fig. 4. — Expected NFW concentrations as a function of redshift for various mass halos from [Duffy et aT] | |2008l l (see Tabic [5|l . Also 
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Fig. 5. — The effect of cosmology. Left: the Millennium simulation (NetoOT) used WMAPl with a higher erg resulting in higher 
concentrations. Right: 3 different simulations with different cosmologies used. 
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Fig. 6. — Comparing several studies, including both simulated and observed NFW c(Af). 



TABLE 3 

NFW £200(^^200) FIT PARAMETERS (2 = 0): C20O ■ 



CO(M200/Mo)-'»±AC200 



Sample 




Cosmology 


CO 


Mo[MQ/h] 


a 


^ logio C200 


DuffyOS 


all 


WMAP5 


5.74 


2 


X 10^2 


0.097 


0.15 


DuffyOS 


relaxed 


WMAP5 


6.67 


2 


XIO12 


0.092 


0.15 


MaccioOS 


all 


WMAP5 


6.12 




IO12 


0.110 


0.130 


MaccioOS 


relaxed 


WMAP5 


6.76 




IOI2 


0.098 


0.105 


MaccioOS 


all 


WMAP3 


5.24 




IOI2 


0.088 


0.132 


MaccioOS 


relaxed 


WMAP3 


5.87 




IOI2 


0.083 


0.109 


MaccioOS 


all 


WMAPl 


7.57 




1012 


0.119 


0.129 


MaccioOS 


relaxed 


WMAPl 


8.26 




1012 


0.104 


0.111 


Neto07 


all 


WMAPl 


4.67 




10" 


0.11 


0.094 


NctoOT 


relaxed 


WMAPl 


5.26 




1014 


0.10 


0.061 



TABLE 4 



NFW 


Cvir{M^ir) FIT PARAMETERS {z = 0) 


Cvir =Co(M 


„»r/Mo)-"±A 


Cvir 


Sample 




Cosmology 


CO 


Mo[Me/h] 


a 


Alogioc„i^ 


DuffyOS 


all 


WMAP5 


7.96 ±0.17 


2 X IO12 


0.091 ± 0.007 






DuffyOS 


relaxed 


WMAP5 


9.23 ±0.15 


2 X IO12 


0.089irol§ 






MaccioOS 


all 


WMAP5 


8.41 


IO12 


0.108 






MaccioOS 


relaxed 


WMAP5 


9.35 


IO12 


0.094 






MaccioOS 


all 


WMAP3 


7.26 


IO12 


0.086 






MaccioOS 


relaxed 


WMAP3 


8.22 


IO12 


0.080 






MaccioOS 


all 


WMAPl 


10.26 


1012 


0.114 






MaccioOS 


relaxed 


WMAPl 


11.25 


1012 


0.099 






Hennawi07 


all 


WMAPl 


12.3 


1.3 X 10" 


0.13 


0.098 


Gcntile07 


all 


WMAP3 


13.6 


lOii 


0.13 




BuUockOl 


all 


WMAPl 


9 


1.3 X 10" 


0.13 


0.14=1 


Comcrford07 


all 


observed 


14.5 ±6.4 


1.3 X 10" 


0.15 ±0.13 


0.15 



al.| ( |2002| footnote 10) claim that the scatter of A loi 
(2OOII was a bit too high and should actually be 0. 



^ Wechsler et 

Bullock et al. ('JoXUT 

with their own measured scatter. 



°giO Cuir = 0.18 reported by 
.14, thus bringing it in lino 



TABLE 5 

NFW c{M, z) FIT PARAMETERS FOR < 2 < 2: 



c = co(M/Mo)--{l + 2)-'^ 



Sample 




Cosmology 


A 


CO 


Mo[M0/h] 


a 


/9 


DuffyOS 


all 


WMAP5 


200 


5.71 ±0.12 


2 X IO12 


0.084 ±0.006 


0.47 ±0.04 


DuffyOS 


relaxed 


WMAP5 


200 


6.71 ±0.12 


2 X IO12 


0.091 ± 0.009 


0.44 ±0.05 


DuffyOS 


all 


WMAP5 


vir 


7.85t»;il 


2 X IO12 


0.081 ± 0.006 


0.71 ±0.04 


DuffyOS 


relaxed 


WMAP5 


vir 


9.23l»;i^ 


2 X IO12 


0.090 ±0.009 


0.69 ± 0.05 



TABLE 
Observed c{M) fit parameters (2 = 



0): c = co(Af/Mo)-"(l + 2)-'3 











Mo 








Mmin 


-iV-Zmax 




Sample 


Analysis 


A, 


CO 


[h-iMo] 


a 


A logic c 


N 


[ft-lAfQ] 


[ft-iAfQ] 


2 


Comerford07'' 


compilation 


vir 


14.5 ±6.4 


1.3 X 1013 


0.15 ±0.13 


0.15 


62 


5 X 1013 


4 X 1015 


0.003 - 0.89 


Buote07 


X-ray 


vir 


9.0 ±0.4 


1014 


0.172 ± 0.026 




39 relaxed 


6 X IO12 


2 X 1015 


0.0033 - 0.2302 


SchmidtAllen07 


X-ray 


vir 


7.55 ± o.go'^ 


8 X 1014 


0.45±0.12<= 




34 relaxed 


2 X 10i4 


4 X IQis 


0.06-0.7 


OkabeOg 


WL 


vir 


4^5-1-3.91 


1014 


0.41 ±0.19 


0.19 


30 


2 X 10i4 


1.5 X IQis 


0.15-0.30 


Johnston07 


stacked WL 


200 


4.1 ± 1.2 


1.3 X 10" 


0.12 ±0.04 




130,000 


5 X IOI2 


5 X 1014 


~0.25 


MandelbaumOS^l 


stacked WL 


54 


4.6 ±0.7 


1014 


0.13 ±0.07 




222,699 


3 X 1013 


6 X 1014 


~0.22 


MandelbaumOS'' 


stacked WL 


200 


2.5 ±0.4 


6 X 10" 


0.13 ±0.07 




222,699 


3 X IO13 


6 X 1014 


~ 0.22 



Note. — SchmidtAllen07 find /3 = 0.71 ± 0.52=, but all others fix p 

^ Includes BuoteO? and SchmidtAllen07 

'^ Quoted uncertainties are 95% rather than 1-cr 

'I MandelbaumOS used H^ir = 200 (not Ac = 200) and assumed Q,m = 

" Converted from previous line 



1 while fitting only (co,o). 



0.27, corresponding to Ac = f!mA„ 



:54. 
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Fig. 7. — Ratio of NFW C200 for relaxed vs. all 2 = halos in 
various simulations. 




2 



;;; DuffyOS;; ;;;; 
(Einasto fits) 


^^ relaxed (z=0) 
^^ relaxed z=0 

relaxed z=l 

relaxed z=2 

NFW 


"**Os_ 





Fig. 8. — For a fixed NFW curve (rs,ps), illustratio n of conver- 
sion of (c, M) between different values of Ac. From Fig. jllj we find 
Cioo ~ 1.37c2oo arid Mioo ~ I.3M200. The two (perfectly overlap- 
ping) black curves have identical {rs,ps)- Dots mark r^, while the 
dashed lines mark the overdensities and r^ir for Ac = 200 and 100. 
Note that for fixed c, r^ and r^^^ vary with M^i^ (black vs. yellow). 
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Fig. 8. — Expected concentrations derived from Einasto profiles 
com pared to NFW profi les. Plotted are relations for relaxed halos 
fromjpuffy et al.l l|2008||. The yellow NFW lines were plotted in 
Fig.l^p-'^-^ 



TABLE 7 

Simulations 




Fig. 9. — Conversion from NFW C200 to c^ir for Ac = 
94, 100, 120, 140, 160, 180, 200. The relations are extremely linear 
for C200 > 2. 



Simulation 


Cosmology 


A/min 


jWniax 


N 


N 


particles 




(see Table sl 


[h-^M@] 


[h-iMo] 


all 


relaxed 


within r„ir 


DuffyOS 


WMAP5 


10" 


lOis 


1,269 


561 


10,000 


MaccioOS 


WMAP5 


lOio 


1015 


9,988 


7,060 


500 


NetoOT (Millennium) 


WMAPl 


IO12 


1015 


53,626 


39,330 


10,000 


HennawiOT 


WMAPl 


10" 


lOis 


878 






BuUockOl 


WMAPl 


10" 


1014 


~5,000 




150-120,000 


Gcntilc07 (NFW96) 


WMAP3 


3 X 10" 


3x 1015 


19 




5,000-10,000 



TABLE 8 

COSMOLOGICAL PARAMETERS 



Author 


WMAP 


nm 


o-g 


BuUockOl 


WMAPl 


0.3 


1.0 


HennawiOT 


WMAPl 


0.3 


0.95 


NetoOT: Millennium 


WMAPl 


0.25 


0.90 


MaccioOS 


WMAPl 


0.268 


0.90 


MaccioOS 


WMAP3 


0.238 


0.75 


MaccioOS, DuffyOS 


WMAP5 


0.258 


0.796 




WMAP7 


0.26 


0.803 
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Fig. 10. — Parameters a and b for conversion of NFW c^ir ~ ac2oo + b. Dashed lines are the relations given in Eqs.|19|fc|20| 
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Fig. 11. — Ratios of r^ir (left) and M^ir (right) to the A^ = 200 values as a function of NFW (c, A^). Note that Tyir/r2oo = Cyir/c2oo 
since we hold Vg (and ps) fixed and we define r^ir = c^irTs and r2oo = C200^s- 
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Fig. 12. — The relati on b etween r„ir and M„i^ is independent of the assumed profile, only dependent on the chosen overdensity, cosmology, 
and halo redshift (Eq. [23|. For a 2 = halo with r^ir = 1 Mpc/h, M^ir = 1.16 X IO^'^Mq/H (Ac/100). 
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Fig. 13. — Qm{z) = 1 — Qf^{z) in a flat universe with a cosmological constant 
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Fig. 14. — Ac{z): The average overdensity above Pcrit within a collapsed halo 
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Fig. 15. — Ayir(z): The average overdensity above pm = f^mPcrit within a collapsed halo 
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Fig. 16. — Power law slopes 7' for 3-D density p oc r""^ . Here 
we compare the NFW, Einasto, and Sersic profiles. For Sersic we 
use the Prugniel-Siniien approximation. A realistic r ange of a and 
n are plotted for the profiles , as found in sim ulations ( [Merritt et al.j 
|2006||Gao"et al., 2008, .Navarro et al.|2010| l. 
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Fig. 19. — NFW power law slopes for projected mass within a 
cylinder M(< ij), projected mass density k(R), and 3-D density 

p{r). 



Fig. 17. — Power law slopes 7' for 3-D density p oc r~'^ . In this 
figure we add the Dehnen-McLaughlin and Stadel-Moore profiles 
with 7 = 7/9 and A = 0.10, respectively. 



steep 5 



NFW 

Dehnen & McLaughlin 7 
Dehnen & McLaughlin 7 
Dehnen & McLaughlin 7 
Dehnen & McLaughlin 7 



flat 0-7t 




Fig. 18. — Power law slopes 7' for 3-D density p oc r~"' . Here we 
compare various v alue s of 7 for Dehnen-McLaughlin as observed in 
simulations (Merritt et al. 20061. The value marked with a * (7/9) 
is used in their 2-parametcr fitting formula. 
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